Astronomy & Astrophysics manuscript no. 25504_finalversion ©ESO 2015 

February 26, 2015 


A hot bubble at the centre of M81 
(Research Note) 

T.V. Ricci, J.E. Steiner, and L. Giansante 

Institute) de Astronomia, Geofisica e Ciencias Atmosfericas, Universidade de Sao Paulo, 05508-900, Sao Paulo, SP, Brasil 
e-mail: tvricci@iag.usp.br, joao.steiner@iag.usp.br, louise.martins@usp.br 

Received <date> / Accepted <date> 


in 


ABSTRACT 


o 

(N 

X 

OJ 

X 

'xj“ 
(N 


< 

o 


Context. Messier 81 has the nearest active nucleus with broad Hu emission. A detailed study of this galaxy’s centre is important for 
understanding the innermost structure of the AGN phenomenon. 

Aims. Our goal is to seek previously undetected structures using additional techniques to reanalyse a data cube obtained with the 
GMOS-IFU installed on the Gemini North telescope (Schnorr Muller et al. 2011). 

Method. We analysed the data cube using techniques of noise reduction, spatial deconvolution, starlight subtraction, PCA tomography, 
and comparison with HST images. 

Results. We identified a hot bubble with T > 43500 K that is associated with strong emission of [N II]/15755A and a high [O 
I]/16300/Hu ratio; the bubble displays a bluish continuum, surrounded by a thin shell of Ha + [N II] emission. We also reinterpret the 
outflow found by Schnorr Muller et al. (2011) showing that the blueshifted cone nearly coincides with the radio jet, as expected. 
Conclusions. We interpret the hot bubble as having been caused by post starburst events that left one or more clusters of young 
stars, similar to the ones found at the centre of the Milky Way, such as the Arches and the IRS 16 clusters. Shocked structures from 
combined young stellar winds or supernova remnants are probably the cause of this hot gas and the low ionization emission. 
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1. Introduction 


2. Observations and data reduction 




j Messier 81 (NGC3031) is an SA(s)ab galaxy, known to host 
a low ionization nuclear emission region (LINER) (Heckman 
t-H 1980) associated with a type 1 active galactic nucleus (AGN) 
(Peimbert & Torres-Peimbert 1981). At a distance of 3.5 Mpc, 
\o this makes it the closest case of a known type 1 AGN. Radio 
r- emission in the form of a compact core, a nuclear jet (Bietenholz 
et al. 2000; Martf-Vidal et al. 2011), and nuclear X-ray emission 
(Young et al. 2007) are associated with this AGN. 

Although photoionization by an AGN is one of the possi- 
f—s ble causes of excitation for LINERs (Ferland & Netzer 1983; 
Halpern & Steiner 1983), shocks (Heckman 1980) and hot old 
H stellar populations (Binette et al. 1994) have also been proposed 
as potential sources. Given its distance, this is an important ob- 
• "Z- ,j ect f°r studying emission line properties with high spatial reso- 
lution. 


In this research note, we reanalyse a data cube of M81 
(Schnorr Muller et al. 2011) obtained with the integral field unit 
(IFU) of the Gemini Multi-Object Spectrograph (GMOS) in¬ 
stalled on the Gemini North telescope. We report the detection 
of hot (T>40000 K) gas, located ~ 0.3 to 0.8 arcsec (5 to 14 pc) 
southward from the nucleus of this galaxy. We also reinterpret 
the spatial orientation of the nuclear gaseous disc and outflow 
reported by Schnorr Muller et al. (2011) 1 . 


1 Allan Schnorr-Miiller was contacted about this new result in M81 
and also about the correction to the spatial orientation of the data cube. 


M81 was observed with the Gemini North telescope on 2006 De¬ 
cember 31 (programme GN-2006B-Q-94). We used the GMOS- 
IFU in two-slit mode. Three exposures of 530s covered the cen¬ 
tral region of M81, with a field of view (FOV) of 7 x 5 arcsec 2 . 
The R400 grating was used, resulting in an observed wavelength 
range of 5600-7000A and spectral resolution R ~ 1800. The raw 
data were reduced with the standard Gemini package under the 
IRAF environment. Bias, flat-fields, calibration lamps, and spec- 
trophotometric standards were obtained for the overall correc¬ 
tion and calibration of the data. Cosmic rays were removed with 
the LACOS algorithm (van Dokkum 2001). Finally, three flux- 
calibrated data cubes were created with a spatial sampling of 
0.05 arcsec. 

In addition to the basic reduction steps, we also adopted 
additional procedures to improve the quality of the data cubes 
(Ricci et al. 2014b; Menezes et al. 2014). First, we corrected 
each data cube for the effects of differential atmospheric refrac¬ 
tion (DAR) using an algorithm developed by us. Secondly, we 
calculated the median of the three data cubes to avoid bad pixels 
and cosmic rays that had not been properly removed by LACOS. 
We then removed high-frequency noise from the spatial dimen¬ 
sions by convolving each image of the median data cube with 
a Butterworth low-pass filter (Gonzales & Woods 2008), with a 
spatial-frequency cut of 0.2Fvy, where F,vk = 0.5 spaxels 1 is 
the Nyquist frequency and n = 2 (see Ricci et al. 2014b for more 
details on the Butterworth filtering process). Instrumental finger¬ 
prints were removed with principal component analysis (PCA) 
tomography (Steiner et al. 2009; Ricci et al. 2014b; Menezes 


Article number, page 1 of 5page.5 




A&A proofs: manuscript no. 25504_fmalversion 


et al. 2014). Finally, we deconvolved each image of the median 
data cube using a Richardson-Lucy algorithm (Richardson 1972; 
Lucy 1974), with ten iterations, assuming a Moffat point spread 
function (PSF) with full width at half maximum (FWHM) of 1.0 
arcsec and /3 = 2.9. The deconvolved PSF is described by a Gaus¬ 
sian function with FWHM = 0.89 arcsec, measured in the image 
of the broad component of the Her emission line. 

3. PCA tomography 

PCA tomography (Steiner et al. 2009) consists of applying PCA 
to data cubes with the aim of extracting useful information and 
removing instrumental noise (see e.g. Ricci et al. 2011; Menezes 
et al. 2013; Ricci et al. 2014b; Menezes et al. 2014). This method 
allows the separation of information originally expressed in a 
system of correlated coordinates into a system of uncorrelated 
coordinates (eigenvectors or eigenspectra). Tomograms are pro¬ 
jections of eigenspectra on the data cubes and reveal where the 
correlations between the wavelengths occur on the spatial di¬ 
mension of the data cubes. 

Schnorr Muller et al. (2011) applied PCA tomography to the 
data cube of M81 in the 6200-6700 A spectral range. They ar¬ 
gued that the bipolar stmeture seen in the tomogram related to 
the third eigenspectra is associated with a gas disc. We agree 
with their interpretation; however, because of an error in the 
spatial orientation of the FOV (the x-axis of these authors cor¬ 
responds to -x in our figures), the position angle (PA.) of the 
gas disc is now -3°±1°, measured as the angle of the line that 
connects the positions of minimum and maximum weights in 
the tomogram. We show the tomogram and the eigenspectrum 
related to the gas disc in Fig. 1. It is worth noticing that, in ad¬ 
dition to the disc-like kinematics associated with the Hcr+tN II] 
lines, one can also see a strong correlation between [O I] lines 
and the blueshifted side of the disc to the south of the nucleus. 
Moreover, a significant bluish continuum is also related to this 
feature. Such correlations have not been seen in other disc-like 
kinematics of early-type galaxies (Ricci et al. 2014b). 

The tomogram and the eigenspectrum 4, shown in Fig. 2, 
also suggest the motion of the gas, but the direction of this kine¬ 
matic feature is almost perpendicular to the gas disc (PA. = - 
103°+ 4°). This is probably related to an outflow, since a nuclear 
radio jet is pointing to the same direction as the component that 
is in blueshift relative to the nucleus (PA. = 65°: Marti-Vidal 
et al. 2011). Schnorr Muller et al. (2011) also proposed that the 
fourth eigenspectrum is related to an outflow; however, in their 
study, because of the incorrect spatial orientation of the FOV, 
the nuclear radio jet pointed to a structure that was in redshift 
relative to the nucleus of M81. 

We also applied PCA tomography to the 5675-6700 A spec¬ 
tral range. The main difference here is that we included in 
the analysis the [N IIJT5755 emission line, which is sensitive 
to the high temperatures of the ionized gas. Again, the first 
and second eigenspectra are related to the galactic bulge and 
to the AGN, respectively. However, the third eigenspectrum, 
shown in Fig. 3, reveals a correlation between the [N II]T5755, 
[O I]/LI6300,6363, and Ha+[N II]Td6548,6583 emission lines. 
These lines are also correlated with the blue region of the contin¬ 
uum. The tomogram related to this eigenspectrum, also shown in 
Fig. 3, reveals a region associated with these emission lines south 
of the nucleus, again coinciding with the blueshifted side of the 
gaseous disc. The interpretation is that this region is composed 
of very hot gas and that the continuum is bluer than the overall 
stellar component along the FOV, which indicates the presence 
of young stellar populations. 
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Fig. 1. Tomogram and eigenspectrum 3 of the PCA tomography applied 
to the 6200-6700 A spectral range. The anti-correlation between the red 
and the blue wings of the [N II] and Ha emission lines associated with 
a bipolar structure seen in the tomogram is a typical signature of gas 
kinematics. In the tomogram, the stmeture in blue (red) corresponds to 
the blue (red) wings of the [NII] and Ha emission lines. In this case, we 
interpreted this kinematics as a gas disc with P.A. = -3°±1°. The white 
circle marks the position of the AGN as given by the image of the broad 
component of Ha, and the white arrow corresponds to the direction of 
a resolved nuclear radio jet with P.A. = 65° (Martl-Vidal et al. 2011). 


4. Emission line properties 

In order to study emission lines emerging from the hot gas de¬ 
tected with PCA tomography, we subtracted the stellar compo¬ 
nents of each spectrum of the data cube by means of the spectral 
synthesis technique. We used the software STARLIGHT (Cid 
Fernandes et al. 2005) with the stellar population base proposed 
by Walcher et al. (2009), which contains 120 stellar populations 
with ages between 3 and 12 Gyr; abundances of [Fe/H] = -0.5, 
-0.25, 0.0, and 0.2; and [a/Fe] = 0.0, 0.2, and 0.4. The result of 
subtracting the solutions of the spectral synthesis from the spec¬ 
tra of the data cube for each spaxel is a gas cube, i.e. a data cube 
containing only the gas component. Since we are interested in 
analysing only the narrow component of the emission lines, we 
removed the broad component of the Ha line from the gas cube 
of M81 in a way similar to that of Ricci et al. (2014a). 

We measured the [O I]46300/Ha and [N II]45755/[N 
II]/l/l6548+6584 line ratios along the FOV. To do so, we adjusted 
the profile of each line with a Gauss-Newton algorithm to fit non- 
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Fig. 2. Tomogram and eigenspectrum 4 of the PCA tomography applied 
to the 6200-6700 A spectral range. In this case, the gas kinematics is 
probably associated with outflows. We note that the bipolar structure is 
almost perpendicular to the gas disc (RA. = -103°± 4°). In addition, the 
nuclear radio jet is in the same direction (P.A. = 65°) as the component 
of the outflow that is in blueshift relative to the nucleus of M81. 

linear functions. For each spaxel, we assumed that the [O I], [N 
II], and Her lines have the same velocity and FWHM. In Fig. 4 we 
show the [O I]/Hcr line ratio map. We also show the gas temper¬ 
ature derived from the [N II]T5755/[N II]TT6548+6584 (Oster- 
brock & Ferland 2006), assuming a low gas density regime. The 
outermost isothermal curve (red line in Fig. 4) corresponds to T e 
= 43500 K. The maximum emission of [N II]5755 is located at 
0.29 ± 0.07 arcsec south of the nucleus, while the [O I]/Ha ratio 
is at 0.7 arcsec, where it reaches > 0.6, and the maximum of the 
temperature seems to be located at 0.8 arcsec. We conclude that 
the bubble comprises a region from 5 to 14 pc from the nucleus 
of M81. Given the limitations of our resolution it is not possi¬ 
ble to assert its real dimension, although the HST Her emission 
(Fig. 4) suggests that the bubble could extend as far as 2 arcsec 
(35 pc). It would be highly desirable to have observations with 
higher spatial resolution, for example with adaptive optics of the 
[Fe II]/L112567, 16435 lines. 

The [N II]T5755 line has a width of 388 ± 9 km s _1 ; we 
do not have an indication of the gas density because the [S H] 
doublet fell in the spectral gap of the CCD detectors. Therefore, 
no mass estimate can be made. 

In Fig. 4 we also show the ratio between the archival images 
of the narrow Htr+[N II]TT6548+6584 filter (namely F658N) 
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Fig. 3. Tomogram and eigenspectrum 3 of the PCA tomography applied 
to the 5675-6700 A spectral range. This range includes the [N II]T5755 
emission line, which is sensitive to the high temperatures of the ionized 
gas. The eigenvector reveals a correlation between all emission lines 
and the blue component of the continuum. 

and of the I band obtained with the Wide Field Camera (WFC) 
of the Advanced Camera for Surveys (ACS) on board the Hubble 
Space Telescope (HST). This ratio shows the strong emission of 
the Hcr+tN II]/LI6548+6584 lines the peak of which is taken to 
be the centre of the FOV, as it presumably indicates the position 
of the AGN. For the GMOS frame, we assume the centre to be 
derived from the broad Her red wing. With this assumption we 
can overlay images of both the HST and Gemini. In addition to 
the AGN, a faint shell of emission nearly encircling the region 
of hot gas is also seen in the HST image. The maximum temper¬ 
ature is located 0.8 arcsec (14 pc) south of the nucleus. 


5. Discussion and conclusions 

In this work, we reanalyse the data cube of M81 studied by 
Schnorr Muller et al. (2011). The authors found that this galaxy 
has a gaseous disc as well as an outflow associated with the 
AGN. We found that their x-axis was inverted from what it 
should have been and, therefore, that the new P.A. of these 
gaseous features are quite distinct. In particular, the radio jet ob¬ 
served by Bietenholz et al. (2000) and Marti-Vidal et al. (2011) 
seems to be associated with the redshifted cone/outflow struc¬ 
ture (ionization cone) detected by Schnorr Muller et al. (2011); 
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Fig. 4. (Har+[N II])/I ratio using images from the ACS aboard the HST. 
We note that the peak of the emission coincides with the position of 
the nucleus. The magenta contours are curves of isotemperatures. The 
outer contour has T t , ~ 43500K and the temperature increases inward. 
The cyan contours are associated with the [O I]/Hcr ratio, with values 
of 0.6, 0.7, and 0.8 from the outside to the inside. The x sign marks the 
maximum [N II]T5755 emission. A very hot region associated with high 
values of [O I]/Hcr is probably related to shock-heated regions since 
photoionization models cannot produce such a high temperature. The 
resolved nuclear radio jet has a P.A. of 65". 



for this to make sense, one has to assume that a hypothetical 
counter-jet does not exist. 

With our new orientation of the FOV, the blueshifted 
cone/outflow is co-aligned with the radio jet, for which we as¬ 
sume a new P.A., given by Marti-Vidal et al. (2011) and shown 
in Fig. 2. The asymmetry of an one-sided radio jet is now natu¬ 
rally explained by the Doppler boosting effect. 

When we analysed the data cube with the [N II]/15755A line 
(not included in the Schnorr Muller et al. 2011 study), we de¬ 
tected an unresolved hot spot located 0.8 arcsec south of the nu¬ 
cleus. We estimate a temperature T e > 43500 K for the gas in 
this region. The line ratio [O I]/Hcr = 0.8 is exceptionally high, 
even for LINERs (Ho et al. 1997); it is also higher than on the 
nucleus of M81, where the ratio is [O I]/Ho' = 0.6. High temper¬ 
atures, associated with high [O I]/Ho', are typical signatures of 
mechanical heating such as occurs in shock waves produced, for 
example in supernova remnants (SNR) or shocked stellar winds 
(Dopita & Sutherland 1996; Osterbrock & Ferland 2006; Allen 
et al. 2008). An alternative would be a shock wave produced by 
a superwind induced by the AGN. This interpretation, however, 
seems to be at odds with the jet orientation. The bubble is lo¬ 
cated perpendicular to the direction of the jet; this suggests that 
the bubble is not associated with outflows from the AGN. 

The structure detected in M81 may be similar to (but on a 
larger scale) the structure known as Arches Cluster (Yusef-Zadeh 
et al. 2002), located about 30 pc from the centre of the Milky 
Way, or to the young stellar cluster associated with IRS 16, even 
closer to the centre of our Galaxy (see Genzel et al. 2010 for 
a review). In particular, the Arches Cluster shows a structure 
of cold gas surrounding it. At the same time, hot intracluster 
gas was detected by its X-ray emission and was interpreted as 
a combined shock of stellar winds from young (1-3 Myr) stars 
(Yusef-Zadeh et al. 2002). In M81, as in the Arches Cluster, the 
hot bubble is also encircled by a filament of cooler gas emis¬ 
sion (Fig. 4). Here, young stellar populations revealed by PCA 
tomography may be associated with a more massive cluster or. 


perhaps, a combination of many clusters. For example, Durre & 
Mould (2014) detected four young stellar clusters encircling the 
AGN of the galaxy NGC 2110. Using IFU data with adaptive 
optics in the IR, they associated [Fe II] emissions with shocks 
caused by winds and/or SNR emerging from these clusters. In 
fact, they showed that the peak of the [Fe II] emission in the cir- 
cumnuclear region of NGC 2110 is not collocated with the stellar 
clusters, but lies between them, which indicates that the shocks 
occur in the intracluster medium. In addition to the shocked gas 
from young stellar winds, shock waves from SNR of somewhat 
older (age < 70 Myr) stellar populations may also be at play in 
M81. 

Since M81 is a nearby galaxy, it is possible to discuss the im¬ 
portance of shocks in LINERs even when an AGN is clearly de¬ 
tected. Because of the proximity of M81, we were able to resolve 
the AGN and a region dominated by shocks. If the same sce¬ 
nario applies to relatively distant galaxies, both regions would be 
unresolved and probably classified as pure LINER-like AGNs, 
i.e. with no evidence of shocks. The case of M81 suggests that 
shocks may have an important role even in LINERs that are 
clearly ionized by other sources, for instance AGNs or post-AGB 
stars. 
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